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Abstract 

Stents are medical devices designed to modify blood flow in aneurysm sacs, in order to 
prevent their rupture. Some of them can be considered as a locally periodic rough bound- 
ary. In order to approximate blood flow in arteries and vessels of the cardio-vascular system 
containing stents, we use multi-scale techniques to construct boundary layers and wall laws. 
Simplifying the flow we turn to consider a 2-dimensional Poisson problem that conserves es- 
sential features related to the rough boundary. Then, we investigate convergence of boundary 
layer approximations and the corresponding wall laws in the case of Neumann type boundary 
conditions at the inlet and outlet parts of the domain. The difficulty comes from the fact that 
correctors, for the boundary layers near the rough surface, may introduce error terms on the 
other portions of the boundary. In order to correct these spurious oscillations, we introduce 
a vertical boundary layer. Trough a careful study of its behavior, we prove rigorously decay 
estimates. We then construct complete boundary layers that respect the macroscopic bound- 
ary conditions. We also derive error estimates in terms of the roughness size e either for the 
full boundary layer approximation and for the corresponding averaged wall law. 

Keywords: wall-laws, rough boundary, Laplace equation, multi-scale modelling, boundary lay- 
ers, error estimates. 
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1 Introduction 



A common therapeutic treatment to prevent rupture of aneurysms, in large arteria or in blood 
vessels in the brain, consists in placing a device inside the aneurysm sac. The device is designed 
to modify the blood flow in this region, so that the blood contained in the sac coagulates and the 
sac can be absorbed into the surrounding tissue. The traditional technique consists in obstructing 
the sac with a long coil. In a more recent procedure, a device called stent, that can be seen as 
a second artery wall, is placed so as to close the inlet of the sac. We are particularly interested 
in stents produced by a company called Cardiatis, which are designed as multi- layer wired struc- 
tures. Clinical tests show surprising bio-compatibility features of these particular devices and one 
of our objectives is to understand how the design of these stents affect their effectiveness. As 
stent thicknesses are small compared to the characteristic dimensions of the flow inside an artery, 
studying their properties is a challenging multi-scale problem. 

In this work we focus on the fluid part and on the effects of the stent rugosity on the fluid 
flow. We simplify the geometry to that of a 2-dimensional box f2 e , that represents a longitudinal 
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cut through an artery: the rough base represents the shape of the wires of the stent (see fig. [3J 
left). We also simplify the flow model and consider a Poisson problem for the axial component 
of the velocity Our objective is to analyze precisely multi-scale approximations of this simplified 
model, in terms of the rugosity. 

In [4] we considered periodic inflow and outflow boundary conditions on the vertical sides 
ri n ur out of fl e . Here, we study the case of more realistic Neumann conditions on these boundaries, 
which are consistent with the modelling of a flow of blood. 

As a zeroth order approximation to u e , we consider the solution u° of the same PDE, posed 
on a smooth domain fl° strictly contained in fl e . We introduce boundary layer correctors j3 and 
t that correct the incompatibilities between the domain and u . These correctors induce in turn 
perturbations on the vertical sides r in U r out of fl e . We therefore consider additional correctors £j n 
and £ ut, that should account for these perturbations (see figQJ. We also introduce a first order 
approximation, defined in fi , that satisfies a mixed boundary condition (called Saffman- Joseph 
wall law) on a fictitious interface r° located inside fl e . 



Macroscopic problems 



Multi-scale approximations 




Wall law 



Zero order approx 



Micro cell pbm 

Z+ U P, [3 



Full boundary layer 



Micro pbm 1/4 plane 



Figure 1: The exact solution, the multi-scale framework and wall laws 



For the case of Navier-Stokes equations and the Poiseuille flow the problem was already consid- 
ered in Jager et al. [TI] [TO] but the authors imposed Dirichlct boundary conditions on r in U r out 
for the vertical velocity and pressure. Their approach provided a localized vertical boundary layer 
in the e-close neighborhood of r; n U r out . A convergence proof for the boundary layer approxima- 
tion and the wall law was given wrt to e, the roughness size. These arguments are specific to the 
case of Poiseuille flow and differ from the general setting given in the homogenization framework 
|19j . In this work, we address the case of Neumann boundary conditions, where the above methods 
do not apply. The difficulty in this case, stems from the 'pollution' on the vertical sides due to 
the bottom boundary layer correctors. 

From our point of view, the originality of this work emanate from the following aspects : 

- the introduction of a general quarter-plane corrector £ that reduces the oscillations of the 
periodic boundary layer approximations on a specific region of interest. Changing the type 
of boundary conditions implies only to change the boundary conditions of the quarter-plane 
corrector on a certain part of the microscopic domain. 

- the analysis of decay properties of this new corrector: indeed we use techniques based on 
weighted Sobolev spaces to derive some of the estimates and we complete this description by 
integral representation and Fragmcn-Lindclof theory in order to derive sharper L°° bounds. 

- we show new estimates based on duality on the traces and provide a weighted correspondence 
between macro and micro features of test functions of certain Sobolev spaces. 

The error between the wall-law and the exact solution is evaluated on 0°, the smooth domain 
above the roughness, in the L 2 (f2°) norm. This relies on very weak estimates [18] that moreover 
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improve a priori estimates by a y/e factor. While this work focuses on the precise description of 
vertical boundary layer correctors in the a priori part, a second article extends our methods to 
the very weak context [16j in order to obtain optimal rates of convergence also for this step. 

The paper is organized as follows : in section^ we present the framework (including notations, 
domains characteristics and the toy PDE model under consideration), in section [3l we give a brief 
summary of what is already available from the periodic context I it that should serve as a basis for 
what follows, in section |4] we present a microscopic vertical boundary layer and its careful analysis 
in terms of decay at infinity, such decay properties will be used in section [5] in the convergence 
proofs for the full boundary layer approximation as well as in the corresponding wall law analysis. 



2 The framework 

In this work, O e denotes the rough domain in M 2 depicted in fig. O f2° denotes the smooth one, T e 
is the rough boundary and T° (resp. T 1 ) the lower (resp. upper) smooth one (see figH]). The rough 
boundary T e is described as a periodic repetition at the microscopic scale of a single boundary cell 
P°. The latter can be parameterized as the graph of a Lipschitz function / : [0, 2vr[— >] — 1 : 0[, the 
boundary is then defined as 



P° = {y e [0,2tt]x] - 1 : 0[ s.t. y 2 = f( Vl )}. 



(1) 



Moreover we suppose that / is bounded and negative definite, i.e. there exists a positive constant 
5 such that 1 — 5 < f(yi) < 5 for all y\ 6 [0, 2tv]. The lower bound of / is arbitrary and it is useful 
only in order to define some weight function see section [4] We assume that the ratio between 
L (the width of 0°) and 2ire (the width of the periodic cell) is always a positive integer. We 



X2 = 1 



x 2 = 




Figure 2: Rough, smooth and cell domains 



consider a simplified setting that avoids theoretical difficulties and non-linear complications of the 
full Navier-Stokes equations. Starting from the Stokes system, we consider a Poisson problem for 
the axial component of the velocity. The axial component of the pressure gradient is assumed to 
reduce to a constant right hand side C. If we set periodic inflow and outflow boundary conditions, 
the simplified formulation reads : find u e such that 

- = C, in Q e , 

u% = o, on rur 1 , (2) 

u# is x\ periodic. 

In section [3] we should give a brief summary of the framework already introduced in [4] . Never- 
theless the main concern of this work is to consider the non periodic setting (see section |4j where 
we should consider an example of a more realistic inlet and outlet boundary conditions. Namely, 
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we look for approximations of the problem: find u e such that 

- Au e = C, in n £ , 
u e = 0, on PUT 1 , 



(3) 



<9n 



= o, on r in u r, 



in ^ 1 out ■ 



In what follows, functions that do depend on y = x/e should be indexed by an e (e.g. 
U e (x,x/e)). 



3 Summary of the results obtained in the periodic case [4j 

3.1 The cell problems 

3.1.1 The first order cell problem 

The rough boundary is periodic at the microscopic scale and this leads to solve the microscopic 
cell problem: find (3 s.t. 

- A/9 = 0, in Z+UTUF, 
= -y 2 , onP°, (4) 
(3 is yi — periodic . 

We define the microscopic average along the fictitious interface V : (3 = J- f Q w (3(yi, 0)dyi. As 
Z + U r U P is unbounded in the yi direction, we define also 

D 1 - 2 = {v e Ll oc (Z + UTUP)/Dv G L 2 (Z + UTUP) 2 , v is y 1 - periodic }, 

then one has the result : 

Theorem 3.1. Suppose that P° is sufficiently smooth (f is Lipschitz) and does not intersect T . 
Let P be a solution of (j4|), then it belongs to D 1,2 . Moreover, there exists a unique periodic solution 
i] G H 2 (r), of the problem 

< Sr/, (i >=< I, [i >, V/j£ (r), 

where <,> is the (H~ 2 (r), H 2 (r)) duality bracket, and S the inverse of the Steklov-Poincare 
operator. One has the correspondence between (3 and the interface solution n : 

13 = H z+V + Hp? 1: 

where H z +i] (resp. Hprj) is the y\-periodic harmonic extension of rj on Z + (resp. P). The 
solution in Z + can be written explicitly as a power series of Fourier coefficients of 77 and reads : 

H z+V = f3(y)= E me lkyi - ]k]v \ Vy G Z+ , m = n{ yi )e- lk ^d yi , 
, Jo 

k— — oo 



In the macroscopic domain fl° this representation formula gives 



(5) 
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3.1.2 The second order cell problem 

The second order error on T e should be corrected thanks to a new cell problem : find 7 G D 1 ' 2 
solving 

- At = 0, Vy G Z + U T U P, 

r = -y 2 2 , Vy 2 eP°, (6) 
t r periodic in y\ . 

Again, the horizontal average is denoted r. In the same way as for the first order cell problem, 
one can obtain a similar result: 

Proposition 1. Let P° be smooth enough and do not intersect T. Then there exists a unique 
solution t of © in B X - 2 {Z + UTUP). 

3.2 Standard averaged wall laws 
3.2.1 A first order approximation 

Using the averaged value (3 defined above, one can construct a first order approximation u 1 defined 
on the smooth interior domain f2° that solves : 

- Au 1 = C, V.t e 0°, 

u l = e~p^-, Vzer , ^ = 0, Vxer 1 , (7) 
0x2 

u 1 is xi — periodic on Ti n U r out , 
whose explicit solution reads : 



K ' 2 V l + e/3 1 + epJ 



(8) 



Under the hypotheses of theorem 13. 1[ one derives error estimates for the first order wall law 

ll u #- ul IL^°) ^ Kel - 

3.2.2 A second order approximation 

In the same way one should derive second order averaged wall law u 2 satisfying the boundary 
value problem : 

- Au 2 = C, Mx G 0°, 

u 2 = 0, Vi £ T 1 , u 2 is x\ — periodic on r; n U r out , 
whose solution exists, is unique [1] and writes : 

2 . , C ( 2 x 2 (l + e 2 r) e(p-er) 



Now, error estimates do not provide second order accuracy, namely we only obtain 

ll u #- u2 IL^) < k ^i 

which essentially comes from the influence of microscopic oscillations that this averaged second 
order approximation neglects. Thanks to estimates (J5J), one sees easily that these oscillations 
account as if not included in the wall law approximation. 
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3.3 Compact form of the full boundary layer ansatz 



Usually in the presentation of wall laws, one first introduces the full boundary layer approximation. 
This approximation is an asymptotic expansion defined on the whole rough domain Q £ . In a further 
step one averages this approximation in the axial direction over a fast horizontal period and derives 
in a second step the corresponding standard wall law. 

Thanks to various considerations already exposed in [3], the authors showed that actually a 
reverse relationship could be defined that expresses the full boundary layer approximations as 
functions of the wall laws. Obviously this works because the wall laws (defined only on are 
explicit and thus easy to extend to the whole domain Indeed we rc-definc 

C e/3 

ul {x) = — {{l-x 2 )x 2 X[no}+X2X[n'\n ]) ~ =(1-2=2), V.t e fi e (11) 

^ 1 + ep 

while we simply extend u 2 using the formula (|10p over the whole domain. This leads to write : 

dv 1 



1,00 1 , 
= u -f 



du 2 , „, / /a;\ — \ e 2 d 2 u 2 



„}- = U » + e ^ (a!l ,0)^^j-/9j + ^( a : 1 ,0) 



(12) 



For these first order and second order full boundary layer approximations one can set the error 
estimates [4]: 

\\u% - ul'°°\\ L2m < Kei, ||<4 - i#°°|| i3(no) < Ke-i. 

Note that the second order full boundary layer approximation is very close to the exact solution 
in the periodic case, an important step that this work was aiming to reach is to show how far this 
can be extended to a more realistic boundary conditions considered in ([3]). Actually, convergence 
rates provided hereafter and in [16j show that only first order accuracy can be achieved trough 
the addition of a vertical boundary layer (see below) . For this reason we study in the rest of this 
paper only the first order full boundary layer and its corresponding wall law. 



4 The non periodic vertical corrector 

The purpose of what follows is to extend above results to the practical case of We should 
show a general method to handle such a problem. It is inspired in a part from the homogenization 
framework already presented in []I)J[T7] for a periodic media in all directions. The approach below 
uses some arguments exposed in [5] for another setting. 

4.1 Microscopic decay estimates 

In what follows we mainly need to correct oscillations of the normal derivative of the first order 
boundary layer corrector (3 on the inlet and outlet r; n ur ou t. For this sake, we define the notations 
II := U^[Z+ UIUP + 2nke 1 ], the vertical boundary will be denoted E := {0}x]/(0), +oo[ and 
the bottom B := {y e P° ± 2knei} (ci. fig EJ . In what follows we should denote II' := R\, 
B' := R+ x {0} and E' := {0} x R + . 

On this domain, we introduce the problem: find £ such that 

!- A£ = 0, in II, 
f£(0,!/a) = fVifc), onE, (13) 
an an 
f = 0, on B. 

We define the standard weighted Sobolev spaces : for any given integers (n,p) and a real a set 
WZ< p (n) := {v g V'{n)/\D x v\(l +p 2 ) Si±1 ^ IL g Lf(fi), < |A| < 77.} 



6 





tj2 



E' 



n' 



B' 



Figure 3: Semi infinite microscopic domains: H, the rough quarter-plane and IV, the smooth one 



where p := \Jy\ + (2/2 + I) 2 - In what follows we should distinguish between properties depending 
on p which is a distance to a point exterior to the domain II and r = vJ/f+y! the distance to 
the interior point (0, 0). These weighted Sobolev spaces are Banach spaces for the norm 



u\\wt^)--=\ E \\( i +p 2 y 

yO<\X\<m 



the semi-norm being 



Lp(Q) 



p 



We refer to [HI [Hill] for detailed study of these spaces. We introduce a specific subspacc 

W% n (It) = {»e WP' n (U) s.t. v = on B} . 

We begin by some important properties satisfied by £ that will be used to prove convergence 
theorems 15. 1[ 15.21 and 15.31 Such estimates will be obtained by a careful study of the weighted 
Sobolev properties of £ as well as its integral representation through a specific Green function. 

Theorem 4.1. Under the hypotheses of theorem \3.1\ there exists a unique solution of problem 
(|13[) . Moreover £ G W^ ,2 (II) with a £] — ao,ao{ where a$ :— (\p2j~K) and 



\Z(v)\ < 



K 



P(y) 1 ' 



— , Vy G R+ s.t. P >1, 



oyi 



dyi < y+2£ > V ^2 G R+- 

y 2 



where M is a positive constant such that M < 1/(1 — 2a) ~ 10. 

The proof follows as a consequence of every result claimed until the end of subsection 14.1 
Lemma 4.1. In W^' (It) the semi-norm is a norm, moreover one has 

1 



lei 



W«> 2 {IV) - 



2ao 



Va G 



On the vertical boundary E one has the continuity of the trace operator 



k e | 



wi' (E) 



the weighted trace norm being defined as 



i 2 



W 2 ' (dlt) = < u G V (dlt) s.t. 



(l + p 2 )i 



r G i 2 (5n), 



an? 



\y- y'\ 2 



ds(y)ds(y') < +00 
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and 

u £ W^ ,2 (dU) (1 + p z )^u £ W ^"(dU). 

The proof is omitted: the homogeneous Dirichlct condition on B allows to establish Poincare 
Wirtingcr estimates ([6], vol. I page 56) in a quarter-plane containing II. Nevertheless, similar 
arguments are also used in the proof of lemma 14.41 

Lemma 4.2. The normal derivative g := J^(0, 2/2) *s a linear form on Vi^' 2 (II) for every a£l. 

Proof. In Z + , the upper part of the cell domain, the harmonic decomposition of (3 allows to 
characterize its normal derivative explicitly on E' . Indeed 

9 = n { 13 ikrjke- ]kly2 j X[E'] + 9- =■ 9+ + 9- 

\k— — 00 ) 

where g- is a function whose support is located in y2 G [/(0),0]. One has for the upper part 



J Et 9ly2dy 2 <K\\r)f H i {r) , Vae 



thus g+ is in the weighted L 2 space for any power of (1 + p 2 ) 5 , it is a linear form on W2' a (E). For 
<?_, we have no explicit formulation. We analyze the problem ^ but restricted to the bounded 
sublayer P. We define /3_ to be harmonic in P satisfying /3_ = n on T, where n is the trace 
on the fictitious interface obtained in theorem 13.11 and /?_ = — y 2 on P°. Note that thanks to 
standard regularity results n £ C 2 (T) because T is strictly included in Z + U T U P, [7J. So (3 solves 
a Dirichlet j/i-periodic problem in P with regular data. As the boundary is Lipschitz j3 £ H 1 {P) 
and so on the compact interface E c := {0} x [/(0), 0], d n /3 is a linear form on H 2 functions. Then 
because S c is compact : 

d(3 



df3 a < 
-^-vdy 2 < 

Ec dn 



On 



v\\„l,„ , < K\\v\\ i.2 , Vael 

□ 



Lemma 4.3. If a £] — ao : ao[ there exists £ £ M^' 2 (II) a unique solution of problem (j 13[> . 
Proof. The weak formulation of problem (|13[) reads 

(V£,W)n = (<?,«) £ , V«eC°°(IT), 

leading to check hypothesis of the abstract inf-sup extension of the Lax-Milgram theorem [18j [2] , 
for 



p p+00 
/ Vu • Vvcfr/, = / 

Jn+ Jf(o) 



df3 



2- 



a(u, u)= / \7u-Vvdy, l(v) = / j—vdy 

/(o) c*n 

By lemma |4~2I ^ is a linear form on W^' (II). It remains to prove the inf-sup like condition on the 
bilinear form a. For this purpose we set v = up 2a and we look for a lower estimate of a(u, v). 



a(u, up 2a 
Using Holder estimates one has 



/ Vu ■ V (up 2a ) dy = \u\ 2 w i,2 (n) + 2a / p^^uVu -V pdy 




p 2a - 1 uVu-Vpdy<[ j p 2a {-) dy) (/ p 2a \Vu\ 2 dy\ < J- |«|^.- (n) 



7 1,2, 



In this way one gets 



a(u,up Q ) > (1 )\ u \wt 2 (n) 



and if a < ao the inf-sup condition is fulfilled, the rest of the proof is standard and left to the 
reader [3J. □ 
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Thanks to the Poincarc inequality in II \ II' with a = 0, we have 
Corollary 4.1. If a function £ belongs to W ' (II) it satisfies £ G L 2 (B'). 

To characterize the weighted behavior of £ on P' we set uj a (yi) = (yf + 1)~ yi and we give 

Lemma 4.4. If £ in W^ 2 {U) then £ G L 2 {B',uj a ) := {u G V'(B') s.t. f B , £, 2 uj a d yi < oo}. 

Proof. II is contained in a set K + x {— l,+oo}. We map the latter with cylindrical coordinates 
(p,0). Every function of W ' (II), extended by zero on the complementary set of II, belongs to 
the space of functions vanishing on the half-line 8 = 0. Using Wirtinger estimates, one has for 
every such a function. 



i 2 f^arcsin ( i ) 







G) 


dB 



pdOdp < K\\£\\l,i,2 



wr(n) 

because on B' , pdp = yidyi, one gets the desired result. □ 

In order to derive local and global L°° estimates we introduce in this part a representation 
formula of £ on II'. As long as we use the representation formula below, x will be the symmetric 
variable to the integration variable y. Until the end of proposition [2] both x and y are microscopic 
variables living in II. 

Lemma 4.5. The solution of problem (|13|) satisfies £(y) < Kp~ 1+ ^r for every y G II' such that 
p(y) — 1- The constant M can be chosen such that M < 1/(1 — 2a) ~ 10. 

Proof. We set the representation formula 

£(x) = f T x g{y 2 )dy 2 + / ^-^ yi ,0)d yi =: N(x) + D(x), Mx G H', (14) 
7s' Jb< on 

where the Green function for the quarter-plane is 

F x (y) = — (In I a: - y\ + In |x* - y\ - In |a;* - y| - In |x - y|) , 

47T 

with x = (xi,x 2 ),x* = {—x\,X2),x* = (xi, — x 2 ), x = (—x\, —x 2 ). 

The Neumann part N(x). We make the change of variables x = (r cos r sin tf) which gives 



m m <>oo 

lim V 77 fe 7V fe = hm V ^ / fce"^ 2 (ln(s? + (2/2 - z 2 ) 2 ) - ln(x 2 + (y 2 + x 2 ) 2 )) dy 2 

lS j.it»jf».-*( h (i-^ + (?)")- h (i 



— +(7) 

where c = cosi?, s = sini?. Now we perform the second change of variables t& = e - ^ 2 and get 

The last rhs is independent of k; one easily estimates it using the change of variables y = — lnt/r, 
indeed: 

1 / / hiA 2 \ . f°°, , , . f°° e~ r y . 1 



In \\l - — J \dt= ^ ln(l + y)e- r y r dy = ^ — < 

Now because the fictitious interface T is strictly included in Z+ U T U P, /3 G H? oc (Z+ U T U P) 
and thus ||??||#i(r : = EfcLi l^fcl 2 ^* 2 < +00, one has for every finite m 



m /OO \ 2 / OO -. \ 2 

£i.^i< Em 8 * 3 

fc=l \fc=l / \fc=l / 

the estimate being uniform wrt m one has that N < C|| 7 7ll J F^i(r)/' r - 
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The Dirichlet part D(x). We have, by the same change of variable as above (x := r(cos(#), sin i9) 
r(c,s) : 



D(x) 



%2 

~2tt 



1 



1 



{Vi -xi) 2 +x\ {yi+xxY + xi, 



< — 

irr 



o l-2cf + (f) 



s f°° 

■»! = -/ - 
^ Jo 1 



1 



(l-c 2 ) + (c-f) 



where we suppose that c < 1. We divide this integral in two parts, we set m > 1 



s f m 1 

D{x) < — / 

t |> (1 -c 2 ) + (c- JQ-) 

For Ji one uses the Lf oc inclusions : 

s K 

h < 



(l- C 2) + ( c _a) 



(/i+/ 2 )(x). 



rlK(->o)|| L i (0 , m) 



< 



while for I2 one uses the weighted norm established in lemma l4~4l 



L 2 (B' 



< 



2s 



< 



2s 




1 



1 

(l-c 2 ) + (c-f)' 
1 



2/1 



L 2 (B',io a ) 



Vi 2a dyi UWl^b',^) 



2M 



(1-C 2 




L' 2 {B',uj a )' 



where M and M' are Holder conjugates. We choose M'2a > 1 such that the weight contribution 
provided by £ is integrable, this implies that M < 1/(1 — 2a) ~ 10. One then recovers easily 



/ 2 < 



2sA' 



^ V(l-c 2 ) 2M " 



2A" 



(7TX2) 1 " 



We could shift the fictitious interface T to T— Se-2 and repeat again the same arguments because the 
rough boundary does not intersect it. Note that in this case we could establish again the explicit 
Fourier representation formula for (3 and its derivative as in theorem l3.ll Thus we can obtain that 
£ < c(x2 + <5) 1_1 /( 2M ) which shows that £ is bounded in II'. So that on E', one has = 
|£|(1 + x 2 ) 1 ~ 2 ~ M < d . Here one applies the Fragmen-Lindelof technique (see [3], lemma 4.3, p. 12). 
We restrict the domain to a sector, defining Us = n n S, where S = {{p,0) e [l,oo] x [0,7r/2]}. 
On Us, we define w := — 1 + 1/(2M) and v := p ro sm(m8) the latter is harmonic definite positive, 
we set w — £/v which solves 

2 

Aw H — Wv ■ Vw = 0, in Sn, 

v 

with w = on Bs = dS n B, whereas w is bounded uniformly on Es := EDS. Because by 
standard regularity arguments £ € C 2 (n), w is also bounded when p — 1. Then by the Hopf 
maximum principle, we have 

sup \ w\ < sup w < K 
n s an s 

which extends to the whole domain n^, the radial decay of £. □ 
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Deriving the representation formula (|14p , one gets for all x strictly included in II' that 

r r Q qy 

d x £{x)= d xl T x gdy 2 + ——±^ yu 0)d Vl =: N Xl (x) + D Xl (x), Vx E W 
Je> Jb> oxx dy 2 

Lemma 4.6. For any x G II', the Neumann part of the normal derivative d Xl t; satisfies N Xl (x) < 
Kr~ 2 for all x in II' 

Proof. On E' the derivative wrt x\ of the Green kernel reads 

1 1 



X\ 



x\ + (x 2 - y 2 ) 2 %x + (x 2 + V2¥ 



thus using the cylindrical coordinates to express x = (r cos$, r sini?) =: r(c, s) for < s < 1 one 
gets 



pOC ( \ 1 \ 



r 



/ ^ rr^ eH%2 ^ as) 

k Jo 4s 2 (l-s 2 )+ (l-2s 2 - (^) 



1 fOO 1 A 

<4E / yze-^dy! < £ |^| 2 < \\v\\ H - Hry 

, £1^2 Jo X l X 2 , k X 1 X 2 1 ' 

k k 

This estimate is not optimal since it is singular near x\ = or x 2 — 0. But it provides useful 
decay estimates inside II'. 

It's easy to check that N X1 is harmonic, N X1 = g on E' , and that it vanishes on B' . Because 
on E' g is bounded, by the maximum principle N Xl is bounded. We divide II' \ 5(0, 1) in three 
angular sectors : 

Si = {(r, 0) s.t. r > 1, € (^)Lo = {0, |, |, |} 

For Si and 53 we define Vi := ±p~ 2 cos(2#), i = 1,2 which is positive definite and harmonic, 
while for 52, we set Vi := p _2 sin(2?9), that shares the same properties. For each sector we define 
Wi = N xi /vi, it solves 

2 

Awi H Vi>i • Vwi = 0, in Si, i = 1, . . . , 3 

Vi 

The estimate (fT5")) shows that on each interior boundary dSi , w is bounded while on E' U B' it is 
bounded by the boundary conditions that N X1 satisfies. By the Hopf maximum principle [7], one 
shows that 

SUpWi < SUp W\ < OO, SUpU>2 < SUp W 2 < OO, SUpU>3 < SUp W3 < oo 

Si i?=f S 2 #=f,#=f S 3 i?=f,i5 = f 

which implies that 2V X1 < Kr~ 2 for every y € n' s := II' n S*. □ 

In order to estimate -D^ the latter term of the derivative, we introduce the lemma inspired by 
proofs of weighted Sobolev imbeddings in [T3J Q3] . 

Lemma 4.7. If £ £ W*' 2 (II) with a € [0, l/2[ i/ien its trace on a horizontal interface satisfies 

Jo Jo d yi dh <U\\wi- 2 (u) ( 16 ) 

The proof follows ideas of theorem 2.4' in [T3] p. 235, we give it for sake of self-containtness. 
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Proof. We make a change of variables x\ = rcos'd^x^ = rsinz?, j? = 0, leading to rewrite / as 



OO />00 



^0 



|e(r + fe,0)-^(r,0)| ; 



note that the second space variable for £ is now d — 0. We insert intermediate terms inside the 
domain, namely 



KK 



|£(r + /i,0)-f(r + Mtan£)| , |£ (r + Mtan£) - f (r,ataa£) 



-dr dft 



OO /"OO 



|g(r,atan£) -g(r,0)| 



-drdh ^ =: Ji+Ja+i* 



Obviously the terms 2i and I[ are treated the same way. We make a change of variable (r, h 
r tand) 



h 



|£(r(l + tani?),0) - f(r(l + tant?),0)| r 



rdrdd. 



(rtantf) 2 ~ 2a 1 + -d 2 

In order to eliminate the dependence on d in the first variable of £, we then make the change of 
variable (r = f/(l + tan d),d) which gives 



h = 



r^-r; TOT ■ 



^2-2a 



o Jo - - tan 2 - 2Q tf(l + tf 2 ) Jo Jo 

We are in the hypotheses of the Hardy inequality (see for instance [13], p. 203 estimate (7)), thus 
we have 



h < 



(1 - 2a) 2 J 



/ 2 i 



Mfdf 



< Ml*™ 



w„"*(n)! 



in the last estimate we used that in II', the distance to the fixed point (0,-1) can be estimated 
as p 2a := (y 2 + (1/2 + 1) 2 ) Q > (yf + y 2 ) a ='■ r 2a , this explains why we need a positive a in the 
hypotheses. In the same manner 



h < 



o Jo 



|£(r(l + taai>),0)-e(r,0)|' 



(r tan d) 



2-2q 



r dr dd, 



r tan 6 



3£ 

— (r + s,d)ds 



r^dr 



di) 



tan 2 " 2Q d' 



< 



/•tan i9 






/o 


Jo 


dr 



(r(l + er),d)r 1+2a dr 



da 



di) 



tan 2 " 2Q i?' 



where we made the change of variables s = rt and applied the generalized Minkowski inequality 
([13|. p. 203 estimate (6)). Now we set f = r(l + t) inside the most interior integral above 



h < 



tan i9 



dt 




da 


(1 + t) 1 + a 


Jo 


dr 



(r, d)r 1+2 "df 



di) 



tan 2 " 2a d 



Now, we have separated the integrals in i and r, the part depending on t is easy to integrate. Thus 
we obtain 



h < / S(d) 



dr 



(r,i?)f 1+2Q dfdd, where S(d) 



1 



tan 2 " 2 " i? 



(1 + tand) a 



1 



Distinguishing whether tani9 is greater or not than 1, it is possible to show that S is uniformly 
bounded wrt d. This gives the desired result. □ 
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Wc follow similar arguments as in the proof of theorem 8.20 p. 144 in [T^] to claim: 
Proposition 2. Set £ a function belonging to W^' 2 (H), and 



/>oo 

/ G(x,y 1 )£(y 1 ,0)dy 1 , Vx Gil' ', where G(x,y{) 
Jo 



D xi (x) : 

then it satisfies for every fixed positive h 



d dr. r 



dxi dy 2 



yeB' 



/ \D Xl (x u h)\ 2 d Xl < 
Jo 



K 



h l+2 a > 

where the constant K is independent on h. 
Proof. We recall that 

xi - ?/i xi+ j/i 



G := -x 2 



,((xi -yi) 2 +x 2 2 ) 2 ((xi + yi) 2 + x 2 2 ) 2 / 
Because J °° G(x, y\)dy\ = for every x € II' we have 

/>oo 

D Xl (x)= G{x, yi )(Z( yi ,0)-Z(xi,0))dyi, VxelL', 



we underline that G(x, •) is evaluated at x € II' while 0) is taken on B' . By Holder estimates 
in 2/1 with p = 2,p' = 2, we have : 

<J q G\y 1 -x 1 \ dy,^ \ yi _ Xll 2- 2a d Vl (17) 
integrating in x\ and using Holder estimates with p = oo,p' = 1 then 

/ 3 := r\D X1 \ 2 d X1 < sup / ^ij/i-xii 2 ^^! r r^^fe^^i 

Thanks to proposition 14.71 wc estimate the last integral in the rhs above 
h<K sup h{x\, x 2 )\\Z,\\ w i,2 m) i where h{xi, x 2 ) ■= / G 2 \yi 

Jr + 

Considering I4 one has 



xi\ 2 - 2a dy 



t t \ <- is f°° x 2 { Vl -x x f- 2a j x 2 (2/i+x 1 ) 2 (2/ 1 -x 1 ) 2 - 2 " -J 

h{x) - K L i(^WTxW dm+ ((x 1+yi ) 2 + x 2 r dyi 



x 2 {y x -x x f- 2a j X l {yi+Xi y-** 

" A/ o JixT^mWTxW W7+mWTxW " : 



Both terms in the last rhs are treated the same, namely 

/oo 4-2a TS 
_f dz < 
.« (z 2 + l) 4 "4+ 2 " 

which ends the proof. □ 

Remark 4.1. 77m is one of the key point estimates of the paper. One could think of using weighted 
properties of £ of lemma \4-4\ instead of the fractional Sobolev norm introduced from proposition 
\4- 1\ in the Holder estimates (|17|) . This implies to transfer the x\-integral on G, then it seems 
impossible to conclude because 

POO POO 

/ / G 2 uj a (y 1 )dy 1 dx 1 = 00, 
Jo Jo 

which is easy to show if one performs the change of variables z\ = 2/1 — £1,2/1 = yi- 
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In this part we study the convergence properties of the normal derivative of £ on vertical 
interfaces far from E. For this sake we call 

Ui = {y £ n, s.t. yi > I}, Ei = { Vl = I, 2/2 S [/(0), +oo[}, B l = {y€B,y l >l}, 

here we redefine the weighted trace spaces of Sobolev type 

(l + cr 2 )3 Jon? \y-VT 



where we define a := \y — (0, /(0))| = \/y\ + (2/2 — /(0)) 2 - Note that the weight is a distance to 
the fixed point (0,/(0)) independent on I. 

1 2 

Proposition 3. Suppose that v G Wq ct (511/) wra£/i v = on Bi then there exists an extension 
denoted R(v) € X>'(II/) s.t. 

\Vm\w<K\M w fr m y 
where K depends only on || /' || ^ . 

Proof. We lift the domain in a first step in order to transform 11/ in a quarter-plane 11/ . 

y = ^Y):=( X2 ^ m) y yefl/:=(M + ) 2 ; 

if we set v{Y2) = v (Y2 + /(0)) = 17(2/2) then the L 2 part of the weighted norm above reads 
2 f 2 /• - 2 

V , f V I V 



Si (1 + a 2 )^ ~ J El (1 + P + (2/2 - /(0)) 2 )5 - 7 {/}xR+ (1 + / 2 + r 2 2 ) 



(I}xR + (1 + O- 2 ) 2 

where <r 2 = Y : 2 + F 2 2 . If v € W" *; 2 (dIIj) and v = on 5/ we know that ([8], p.43 theorem 1.5.2.3) 

W 2 )i 2 ^<+^, 

which authorizes us to extend v by zero on Ei := {Yi = (} x 1, this extension still belongs to 
- 2 - 

Wq 2 ^. (£/). Arguments above allow obviously to write for every v vanishing on B and v defined 
above 

Ik' 1 1 1 2 = IK'll 1 2 . • 

Here wc use trace theorems II. 1 and II. 2 of Hanouzet [3], they follow exactly the same in our 
case except that the weight is not a distance to a point on the boundary (as in [5]) but it is a 
distance to a point exterior to the domain. So in order to define an extension ([3J p. 249), we set 



V(Y) = [ JC(s)v(Y lS + Y 2 )ds, seR, VF e ft*, 

J\s\<l 

*(y) = $ 



'|s|<l 

Yi-l 



(1 + Y? + P)*J' 
where $ is a cut-off function such that 

Suppfce [0: l/4[, $(0) = 1, $G C°°([0,l/4[), 
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and JC is a regularizing kernel i.e. K. G Cq°(] — 1 : 1[) and LlC(s)ds = 1. Then the extension in 
the quarter-plane domain reads 

w(y) = (w)(r l7 r 2 ) - (w)(yi,-r 2 ), y e u h 

which allows to have w(Y\,0) = for all Y\ G R + . According to theorems II. 1 and II. 2 in [5], one 
then gets 

it) 



(l + a 2 )5 

Turning back to our starting domain II; , we set 

R(v) = w{p^ 1 {y)) = w(y 1 ,y 2 - f(yi)), Vy G U t . 
We focus on the properties of the gradient 

1 /' 



/ {AV y R(v),V y R(v))dy = [ |Vyw| 2 dF, where A = 
Jn, Jfi, 

but the eigenvalues of A are 



/' i+cn 



/\2 i 



2 + (/') 2 ± ^2 + 
A± = 2 ' 

the lowest eigenvalue is positive and tends to zero as |/'| increases. The boundary is Lipschitz 
so that 1 1 f'\ | is bounded. Thus there exists a minimum value of A_. All this guarantees the 
existence of a constant <5'(|i /'II oo) > ® sucn tnat 



6' [ \V y R(v)\ 2 dy<K\\v\\ i, 



which ends the proof. □ 

Thanks to the existence of a lift R(v) , we are able to estimate a sort of weak weighted Sobolev 
norm for the normal derivative on vertical interfaces located at y\ = L/e. 

1 2 

Proposition 4. If v £ (9ITi) and v = on Bl then one has 



(an. 



Proof. The function w given in the hypotheses belongs to the adequate spaces in order to apply 
proposition [3j thus there exists a lift R{v) G W 1,2 (IIt) s.t. R(v) — v on 8Hl. Because £ is 

harmonic and belongs to W 1,2 (II), for any ip £ 2?(IIl) and (p\B L = , one writes the variational 

form: 

" 72 



/ V£ ■ Vtpdy = / —pda(y) = / -^-<pdy : 
Jar, Jdn T , dn JE r , dn 



Then by density and continuity arguments one extends this formula to every test functions in 
p G W ' a (HL_) such that ip = on Bl. As the specific lift R(v) belongs to this space one has 



lz, 2 (ru: 



/ ^fvdy 2 = [ V&R(v)dy< (sup-^ / |V£|V Q dj/ ) ||Vi?(u) 

<^||a w i, a(n ^)l|ViJ(i;)|| La( n^)<^ll«ll Wo ^ (ani) 
which ends the proof. □ 
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4.2 Test functions: from macro to micro and vice-versa 

We suppose that v £ H^^) := {u £ H 1 ^)^ = on T e U T 1 } then j(v) £ H?(dVL e ) which 
implies that v £ H^(Ti n U r out ) and that for any corner 

<■& 

H<t))? T < oo, 
o 1 

where x(t) £ T m U r out is a mapping of the neighborhood of the corners. To the trace of v on Tin 
or r ou t, w e associate a trace of a function defined on OILl which is zero on Bl s.t. 



— ,V2 J := v(0,ey 2 ) = v(0,x 2 ), V.t 2 £ [e/(0), 1] and v ( — , y 2 J := 0, y 2 > -, 



then one has the following connexion between the macroscopic trace norm and the microscopic 
weighted one. 

Proposition 5. Under the hypotheses above on functions v and v, 

\\v\\ i , ~ \\v\\ i , , 

11 M i?7(r, a ur-uri) 11 "w^aru)' 

if v = on r e U r 1 (resp. v = on Bl). 

Proof. Thanks to the change of variables x 2 = ey 2 we have that 

/ v 2 (0,x 2 )dx 2 = e / v 2 ( — ,y 2 ) dy 2 < Kesup(l + a 2 )i [ - — „ , rfy 2 , 

7e/(0) 7/(0) V e / El JSi (1 + O" 2 ) 2 



Conversely 



< Ke-\\v\\ 2 1 2 



v 2 h v 2 1 

= / ^7rdy 2 < sup — - / i l dy 2 



El (1+(T 2 )5 7/(0) (1 + CT 2 )2 Bi (1 + CT 2 )2 7/(0) 



< ^ £ ll«ll^(/(o),i) = ^IMli^iv)- 



For the semi-norm the same change of variable provides an equality due to the homogeneity in e 
i.e. 

II 2 f f Kgg) ~ V ( X ^)\ 2 j j / i~|2 

mi, — II i r^; dx 2 dx 9 = \v\ 1 2 

1 W*) 7 7r Lt |x 2 -4| 2 2 M <; 2 (^) 

□ 

Remark 4.2. We insist on the fact that one can associate traces of v either from T[ n or T out to 
v, the weight that one gains in the microscopic norm comes from the scaling from macro to micro 
and not from the vertical position of the macroscopic interface wrt the origin of the domain Cl e . 

5 A new proof of convergence for standard averaged wall 
laws 

5.1 The full first order boundary layer approximation: error estimates 

The periodic boundary layer approximations given in (fT2"|) introduce some microscopic oscillations 
on the inlet and outlet boundaries T m U r out . We define a new full boundary layer approximation 

ul' ca = u 1 + e^-(x 1 ,0)((3-]3-i in -Ut) (^) (18) 
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where we define 



- = u - . Ut - = e 



~ ( X\ — L X 2 



and £ is the solution of problem (|13p . and £ solves the symmetric problem for r out : 
f -Af = 0, inn_ :=U^ 1 {Z + uruP-2 7 rfcei} 

I dp 

= 7T' on E 

On On 

|=0, on B_ := U^L^P - 27rfcei} 
Every result shown for £ in sections above holds equally for £. One easily checks that 



gfrn 

cm 



i /" I- x 2 

e cm 



j <9£out i 

and "5 \T h 

On 



IdS, f L x 2 
e dn 



We estimate the error of this new boundary layer approximation. We denote r\ ,cc 
it solves 



-Arl'°° =Cx[n.\no], on^ 
dn dn \ e ' e 



5r ^ 



on 



9£ /L x 2 \ 
— -, — 1 on Ii, 
on \ e e 



1,00 



(19) 



eg- (xx , 0) (/? - - £ in - Ut) (y'7) =:6 (-'7 



11 1 
e e 



on rVe'°° = 0, on r e 



As is only a first order approximation, a second order error remains in Q e \ fl° . This explains 
the constant source term on the rhs of the first equation in the system above. We then have 

Theorem 5.1. Under the hypotheses of theorem \3.1\ r\ ,oc satisfies 



1,00 1 



\H 1 {i1') 



< e 



Proof. We separate various sources of errors, we set r 1 the solution of the Neumann part of the 
errors, it solves : 



- An = 0, in n e , 

dn _ d£ ( L x 2 
dn dn \ e ' e 
.7-1=0 onPUT 1 , 



on r Q 



dr\ d£ ( L x 2 



dn dn \ e e 



on ry 



(20) 



then the rest r 2 satisfies 



- Ar 2 = C\[n'\n«], in 

— — = on r in u r out , 

an 

r 2 = b — , - 
V e e 



(21) 



on r x ,r 2 = on T £ 



which is the Dirichlct part of the errors and should be evaluated in a second step thanks to 
appropriate extensions and lifts. 
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The Neumann part The variational form of the problem (fSTTj) reads 

/ Vn-Vvdx= f ^(v)dx 2 , VveH 1 ^), 
J si' Jr out dn 

dividing this expression by ||Vu|| i2 (Q«) we first obtain the equivalence : 



sup 

veit^si') 



iVvl 



= Vn 



L 2 (Q= 



L 2 (Sl' 



Indeed, by Cauchy- Schwartz one has easily that the L 2 norm is greater than the suprcmum while 
a specific choice of v = n gives the reverse estimate. Thanks to this, one has 



|Vri|| i2(ne) = sup 



\Vv\ 



L 2 (Sl" 



We underline that we kept the properties of the traces of i/^,(f2 e ) functions inside the sup that 
we aim to evaluate. This norm is lower that the simple -ff~ 2 (r ou t) which authorizes different 
behaviors of test functions near the corners of r out ([5], p 43, thm 1.5.2.3). 
Now the integral in the rhs of the last expression reads in fact : 



r „ 



dri 
dn 



-y(v)dx 2 



r ou1 



d£ ( L x 2 
dn 



r y(v)dx 2 



HO) 



d£ ( L 
dn 



,2/2 l(v)dy 2 , 



where we constructed v as in section I4.2I i.e. v has the same trace as v but v is expressed as a 
microscopic trace function. Thanks to proposition[5]the corresponding microscopic trace v belongs 



i 2, 



to Wq 2 ct (OHl ) and propositions 0] and O give that 
0f ( L 



HO) 



dn 



y 2 j(v)dy 2 < e a \\v\\ T i 



C (9ni) 



<^>IU (ans) <*V>|| ffl( ^ 



The same analysis and convergence rates hold on r; n with the normal derivative of £ ut- All 
together one obtains 

||Vr 1 || La(n . ) <e 1 + a . 

The Dirichlet part We should lift 6, the non homogeneous Dirichlet boundary conditions on 
T 1 defined in (fTi?)) , so we set 

C := b (^-, a;|x [a o], and f 2 := r 2 - <r. 
Standard a priori estimates give 

||W 2 || L 2 (n£) < ||V<r|| L2(ne) +e, 

where the e in the last rhs comes when estimating the constant source term C localized in £ \ f2°, 
indeed : 

{C,v) n *\op < \\C\\ L 2 {n , xna) \\v\\ L 2 {ntXna} < Ve||C|| L 2 (aeVO o ) ||Vw|| i 2 (0£ \ fi o ) < eC'||Vi;|| I ,2 (ne) 
thanks to a Poincarc inequality in the sub-layer. Hereafter we estimate the gradient of the lift, 



|Vc|| L 2 (0e) <eK 



(->-) 



K 



L 2 (r! 



+ e 6n £2 (r 


*) 




L 2 (ri) + 


e Cout L 2 (r!) 


+ e||9 a:i $out|| L 2 (r i 


< Kei + 2 


e 




+ 

L»(ri) 


ae /. n 

9yi \ e ' e / 


L 2 (ri)_ 





< if 



e 2 -2l7 + e !+ Q 



< if e 1+Q 
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where we used the second estimate of theorem 14.11 describing the decay properties of £. This 
ends the proof: the main error is still made when linearizing the Poiscuillc profile in il e \ f2°. As 



„l,oo ._ 



ri + T2 one gets the desired estimate. 



□ 



Here comes the error estimate in the L 2 norm that add approximately an y/e factor to the a 
priori estimates above. 

Theorem 5.2. Under the hypotheses of theorem \3.1[ r\ ,oc , the solution of system (|19p satisfies 



„i>°° I 



Il 2 (o<> 



< Ke 



l+a 



where the constant K is independent on e and a < y/2/n ~ 0.45. 

Proof. We define v a regular solution on the "smooth" domain f2° (in the sense: not rough, in 
particular, f2° is a rectangle) of the problem 



Av = F, 



dv , , 

— - o, on iy u r; ut 
v = o, on r° u r 1 

where the function F is in L 2 (f2°). Thanks to [5] theorem 4.3.1.4 p. 198, one has that there are 
no singularities near the corners i.e. v G H 2 (fl°) n iJ^(ri°) and 



\ v \\H 2 {n«) - ^ll^llL 2 (n )- 



Testing r^'°° against F, one gets 



(rl'°°,F) 



n" 



drl 



dn 



,l,oo 



dv 
' dn 



r!ur e 



where the brackets stand for the duality pairing between H 1 (r in ' U r^ ut ) and -ffg (]?;,/ U r„ ut ), 
whereas the left bracket denotes the standard L 2 (T° U T 1 ) scalar product. This leads to write : 



li 2 (O0) 



< 



dr. 



l,oo 



dn 



H-Mr.n'ur^) 



| i 2( r ur i 



the latter term of the rhs is classically estimated through Poincare on the sublayer and the a priori 
estimates above for the T part : 



li 2 (r 



< e 2 . 



On r 1 there is an exponentially small contribution of the periodic boundary layer and an almost 
e 2 term coming from the vertical correctors £j n and £ out : 



-i>°° I 



lL 2 (ri) 



< Ke 



■ 1 

— ; 

e e 



< 2M 



L 2 (0,L) 



I\ follows using the same arguments as in the proof of the a priori estimates. The astuteness 
resides in the fact that 



h < 



sup 

v£H? (r ou 



< Ke 1+a 



[ h 2 (r out ) 



Indeed HQ(T' out ) functions when extended by zero on r out are a particular subset of -ff^(T out ) 
functions vanishing on <9r out . At this point one uses the same estimates as in the previous proof 
to obtain the last term in the rhs. 

□ 
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5.2 The standard averaged wall law: new error estimates 

We use the full boundary layer approximation above as an intermediate step to prove error esti- 
mates for the wall law. We denote r\ := u e — u 1 . 

Theorem 5.3. Under the hypotheses of theorem ] 3. 11 one has 

llr 1 !! <e 1+a 

II' <Ul£3(«°) - 

Proof. We insert the full boundary layer approximation between u e and it 1 



r\ := u'-u 1 =u e - u]' 00 + u\'°° - u 1 



We evaluate the L 2 (Q") norm of I\ 



h<K\e 



\eJ L 2 (n») \eJ L 2 (n°)J 



< K I 6? + e 



L 2 (Q°) 



while the first term is classical and the estimate comes from ([5|), for what concerns the second 
term, we use the L°° estimates of theorem 14.11 and get 



e( x - 

a \e 



dx = e 



o Jo p 



< e sup p~ 
|o,ilxfo, 



0,i]x[0,i] P 



—dy 

M 

dy < Ke 2 -^- s " 



2+5 



where 8" is a positive constant as small as desired. This ends the proof. 



□ 



6 Conclusion 

In this work we established error estimates for a new boundary layer approximation and for the 
standard wall law with respect to the exact solution of a rough problem set with non periodic 
lateral boundary conditions. The final order of approximation is of e 1+a where 1 + a ~ 1.45 which 
is compatible and comparable to results obtained in the periodic case (see [SJ 0] and references 
there in). 

Establishing estimates in the spirit of very weak solution [18] but in the weighted context 
improves the L 2 (fl°) estimates but requires an extra amount of work not presented here. This is 
done in [16) . we perform also a numerical validation illustrating the accuracy of our theoretical 
results. 
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